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to time-dependent inhomogeneous magnetic fields are investigated, and possible experimental ob- 
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Spin of mobile carriers (electrons and holes) playSj. 
active central role in the current spintronics effortsE. 
where electronic properties are determined, to a great 
degree, by the presence of nonequilibrium spin. The spin 
response to temporal and spatial changes of the mag- 
netic environment determines various properties ofj 
devices. In the recently proposed magnetic diodeH'E 
for example, an inhomogeneous magnetic environment 
results from inhomogeneous magnetic doping. How a 
magnetic diode adjusts to the switching of an external 
magnetic field and to changes in the applied bias de- 
pends on the way the mobile carriers reach equilibrium. 
Since inhomogeneous magnetic fields are ubiquitous in 
spintronic devices (mostly due to the presence of mag- 
netic/nonmagnetic interfaces), it is important to under- 
stand nonequilibrium spin evolution in such fields. In 
this paper we investigate in detail the transient behavior 
of conduction electron spins, within a Boltzmann equa- 
tion model. A unique feature of the model is that it is 
exactly soluble, allowing a detailed study of the transi- 
tion from the ballistic to the diffusive regime. We show 
that in inhomogeneous magnetic fields a spatial separa- 
tion between spins, an analogue of the Stern-Gerlach ef- 
fect, occurs before spin relaxation begins, but spin cur- 
rent vanishes much sooner, at times of the order of transit 
times. 

The model we consider is a degenerate electron gas 
(a metal or semiconductor) in a magnetic field with the 
largest component in the z direction, and with a gradient 
in that direction. The field has also transverse compo- 
nents (as required by V • B = 0), which are essential 
in rendering SG with electron beams uselessQ, but which 
do not hinder an effective spin separation of conduction 
electrons (see below). We show that an effective spa- 
tial spin separation, along with a flow of spin, is possible 
within ballistic and diffusive dynamics, demonstrating a 
Stern-Gerlach-like (SG) effect with conduction electrons. 
The formalism we use, linear response theory within the 
Boltzmann equation, has begu. applied earlier in various 
forms in transport in generalErcl, and more specifically for 
spin transport in the framework of conduction electron 
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spin resonanceOO and giant magnetoresistancetiJ. Here 
we apply this formalism to a special case of spin dynamics 
in an inhomogeneous magnetic field, and solve it exactly 
for specific boundary conditions. 

In the presence of a uniform electric field E and in- 
homogeneous magnetic field B, semiclassical dynamics 
of electrons in (nonmagnetic) solids is governed by the 
Boltzmann equation 
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where /kA ^ /kA(r, i) is the distribution function of 
electrons with lattice momentum k (band index is sup- 
pressed) and spin A (1 or f for up and 1 or | for down), 
at point r and time t. The notation for the drift field is 
simplified as Ea = Ez + A(/is/e)9-Bz/9z, where /is is 
the Bohr magneton and the electron g-factor is taken to 
be 2; the fields are oriented in the z direction. Band ve- 
locity Vk = de-^/dfik, with Sk standing for band energy 
(we consider systems with inversion symmetry where 
band energy is spin independent). Two momentum re- 
laxation processes are distinguished in Eq. |^. First, 
spin-conserving momentum scattering with rate l/r (r 
is momentum relaxation time), and leading to a quasi- 
equilibrium distribution /kA (i5/kA = h\-fk\), in which 
spin up and down electrons have different chemical poten- 
tials. Second, spin-flip momentum relaxation with rate 
1/Ti (Ti is spin-relaxation time), and leading to com- 
plete (momentum and spin) equilibrium at the local and 
instantaneous B-field: f^.^ = /o[ek + ^IiBB{z,t)], where 
fo{e) = I/[exp(£ — /i)/A;BT-|- 1] is the Fermi-Dirac distri- 
bution function with chemical potential /i, temperature 
T, and Boltzmann constant ks {Sf^x ~ /kA — /kA)- 
writing the Boltzmann equation as Eq. ^ we neglect 
the Lorentz force as unimportant, as the largest part of 
the magnetic field is oriented along the same direction 
as the drift velocity itself (see below for the reasons why 
we also neglect the orbital effects of the small transverse 
magnetic fields). The dynamics of the transverse spin 
components {x and y) is also not considered, as it is 
masked by their fast precession about B. Finally, the 
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contribution of the electronic magnetization to the mag- 
netic field is neglected for our nonmagnetic systems. The 
relaxation time approximation used in Eq. |^ is good for 
all practical purposes, but some caution is needed es- 
pecially at low temperatures, as shown in Appendix 
Generalization of Eq. |^ to k-dependent r, Ti (which can 
vary more wildly over the Fermi surface than r-see Ref. 
[l^ , and ^-factor is straightforward. 

We search for the solution of Eq. ^ in the form /kA = 
/kA ~ {9fo/d£u)4'k\j and write the quasi-cquilibrium 
distribution function as /kA = /kA ^ {dfo/dei^}fix, 
where the nonequilibrium chemical potential jix = {4>k\) 
is obtained self-consistently from the integral relation 
for spin conservation, X^k'^/kA = 0. The angular 
brackets introduce the Fermi-surface averaging: (...) = 
J2k..-dfo/dek/J2kdfo/dek. After linearizing Eq. |l| in 
terms of 6 and E\ , we obtain 
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Particle number conservation requires that (i 
vanishes. The total spin density is 

S = -gpfJ-sB + (l/2)gF(0kT - '/'ki) 



OkT 



), (2) 



(3) 



»ki; 



(4) 



where gp — — 2^^9/o/9ek is the density of states, per 
unit volume, at the Fermi level. The first term on the 
RHS of Eq. y is the equilibrium spin value, yielding the 
electron gas paramagnetic susceptibility of ^i%gF, while 
the second part, JS*, represents the nonequilibrium con- 
tribution to spin density. The spin current density is 



Js = (l/2).gF(fzk(0kT - 0ki)), 



(5) 



and is connected to S via the continuity equation derived 
from Eq. p|. 



dS/dt -t- dJs/dz = -5S/Ti, 
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which., together with the linear response equation (see 
Ref.tj for a systematic treatment of linear spin trans- 
port). 



J, = -DdSS/dz, 



(7) 



where D = (wfk)''' is the electron diffusivity constant, 
gives the diffusion formula for investigating diffusive spin 
transport. 



d6S d^6S _ SS_ dSo 
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We first study transient phenomena that describe evo- 
lution of 5" and Js towards equilibrium. Consider an 
unpolarized sample (whose band structure is assumed 



isotropic, Vk = h]s./m, where m is electron band mass), 
stretching from to L along the z-axis, with no charge 
current {E = 0). At t — 0, magnetic field B{z) — 
Bn + Biz is apphed. Our goal is to find, by solving Eq. ||, 
{(j)) EE {(j)^) = -{(pi) [so that S = -gpf^BB + gpitj))], sub- 
ject to the initial condition (0(z,O)} = fj.B[B{z)], where 
[i?(z)] — Bo + Bi[z] is the even periodic extension of B{z) 
from interval (0, L) to the whole z-axis; thus formulated 
initial condition guarantees that spin current vanishes at 
the boundary-the assumption well justified in nonmag- 
netic interfaces with negligible spin-flip scattering. The 
spin profile can be written as S = Sh + Sim where the ho- 
mogeneous Sh and inhomogeneous Sm spin components 
are 

Shit) = -gFMBo + lBiL)[l-K{0,t)], (9) 
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with g„ = (2n-|- / L. The sum comes from the Fourier 
expansion of [-6(2)]. Kernel K{q,t) describes the time 
evolution of the Fourier g-components of the nonequilib- 
rium spin: 

K{q, t) = SS{q, t)/SS{q, 0) = (</)(g, t))/{<l>{q, 0)). (11) 

Having the spin, the spin current can be calculated from 
the continuity equation as 
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Equation ^ gives also an exact quasilocal relationship be- 
tween spin current and spin, valid at all times, and ex- 
pressed in terms of the Fourier coefficients as 



Js (g, t) = ivpSSiq, 0)i?i (q, t) + ivF - [5S * Ri] (g, i)-(13) 

T 



Here vp \s the Fermi velocity, 
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with X = qvpt and 1/t„i = 1/t + 1/Ti the total momen- 
tum scattering rate, and the star symbol denotes tempo- 
ral convolution: 



[/l*/2](g,i)= / dt'h{q,t-t')h{q,t'). 
Jo 
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In real space Eq. 13 expresses Jg in terms of deriva- 



tives (in principle of all odd orders-that is why the term 
quasilocal) of SS. In the diffusive regime, at i ^ r, the 
memory of the initial condition is lost, and Eq. |l^ re- 
duces to Eq. ^ An exact generalization of the diffusion 
equation js] is obtained by substituting Jg from Eq. |l^ to 
the continuity equation. The result is 



K{q, t) = i?o(g, t) + -[K* Ro] {q, t), 

T 



(16) 



where now 

i?o(g,0 = e:>ip{-t/Tm)sm{qvFt)/{qvFt). (17) 

At t > r, Eq. |l^ is equivalent to Eq. ||. 

Equation |^ can be solved exactly with the help of 
Laplace transform. The solution is provided in Appendix 
H. The result is 
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where the vertical bars denote the singular (principal) 
part of the Laurent series in terms of qvpT, of the expres- 
sion inside. An alternative formulation for the Kernel is 




K{q,t) = -{T/t)e- 
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where functions En (x) are described in Appendix |^. 

Let us consider the limiting behavior of K{q, t) for the 
ballistic and diffusive regimes. For ballistic transport, 
t T, the evolution kernel, Eq. |lj, reduces to 
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This is the solution of Eq. in the absence of scattering. 
A finite S{z,t) in the baUistic case is solely due to SG 
effect of semiclassical separation of spins. At the left 
boundary, 



<S'baii(0) = -gFlJ.BBiL{vFt/2L), 



(21) 



in the middle 5baii(i/2) = 0, and at the right boundary 
Shai\{L) = — ^baiilO). The spin separation S{L) — S'(O) 
grows linearly with time, reaching its maximum of about 
—gpfJ-sBii at < = r. For the diffusive transport, t ^ t, 
Eq. gives K w A'diff , where 



(22) 



which is also a solution of Eq. [g. A new time scale, 
tr = L^/Dtt^, appears, for the transit time of a diffusing 
electron crossing the sample (tt^^t is called the Thouless 
time). We consider L smaller than the spia diffusion 
length (which can be as large as a millimeterfiJ) , so that 
It <Ti. For t < It the spin density grows diffusively, 
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(23) 



'S'diff(-^-) = -S'diff(0), while 5diff(£/2) = 0; a large spin 
current flows in the middle of the sample. While the spin 
current vanishes at greater times, t > It, (when drift is 
being balanced by diffusion) , an effective spin separation 
remains almost stationary until t = Ti, when spin re- 
laxation establishes equilibrium. Note that the homoge- 
neous component of spin, Sh, evolves towards equilibrium 
with spin-flip processes only, since K{0,t) = exp{—t/Ti) 
at all times. 



FIG. 1: Calculated normalized spin density 5* (top) and spin 
current Js (bottom) for a model sample defined in the text. 
The curves represent profiles at times t = lO"'^*, ...,10~^ s 
(increasing by a decade), and are denoted by a corresponding 
number 1,...,8 (except for a few cases at the top, where the 
trend is clear). The dotted lines represent the initial ballistic 



transport at t = 10 
the longest times. 



s, while the long-dashed lines are for 



An example of a transient evolution of spin and spin 
current is shown in Figs. |^ and |^. We take a model sam- 
ple of size L — 1/im, with realistic electronic parameters 
r = 0.1 ps, D = 0.01 m^s-^ vp = {SD/rf-^ « 5.5 x 10^ 
ms-\ It = 1?IDt{^ « 10 ps, and Ti 10 ns. Mag- 
netic field is normalized to Bq ~ B\L. We evaluate our 
exact solution, Eq. |l^, numerically to obtain spin, and 
then calculate the spin current from the continuity equa- 
tion. The physics that emerges from our calculation, and 
which can be seen on the model example in Figs. ^ and 
^, is the following. There are four time scales to consider 
(Fig. H). (i) In the ballistic regime (i < r), electron spin 
density at the edges begins to grow as ^ as electrons 
with one spin direction after bouncing off the boundary 
decelerate and stay close, while the electrons with the 
opposite spin accelerate in the other direction. Spin cur- 
rent, which is always largest in the middle of the sample, 
rapidly increases to reach its maximum value at t = t 
(see Fig. ||). Note that positive Js means negative Js, 
and largely a drift spin flow (spin diffusion acts in the 
other direction), (ii) The diffusive regime {t < t < tr) 
is characterized by a further build-up of spin density at 
the edges of the sample (by diffusion) at the rate ~ 
This is accompanied by a decay of spin current, as the 
initial drift is now being balanced by diffusion. We call 
this diffusive SG effect, (iii) The quasiequilibrium regime 
(tx < t < Ti), where momenta are in equilibrium char- 
acterized by a finite difference in the chemical potentials 
of spin up and down electrons, i5S' is spatially uniform 
so that no spin currents flow, and spin densities remain 
almost constant in time, (iv) Finally, in the spin relax- 
ation regime {t > Ti), the uniform nonequilibrium spin 
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FIG. 2: Calculated time evolution of normalized spin density 
S at z = (a), z = L/2 (b) and z = L (c). The dashed line is 
the time evolution of normalized spin current Js at z = L/2. 
The three vertical lines separate the ballistic {t < t — 10~^''s), 
diffusive {t < t < tr ~ lO"^'^ s), quasiequilibrium [tr < t < 
Ti = 10~* s), and spin relaxational {t > Ti) regimes. 
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FIG. 3: Calculated dimensionless spectral response Ri(ijj) + 
R2{z,ui) at z = for the model sample. The two shoulders 
in the real part of the response (and the corresponding peaks 
in the imaginary part) correspond to spin relaxation {lo = 
1/Ti = 10* s"^) and diffusion {uj = l/fr ~ lO" s'^). 



density vanishes and complete equilibrium is established. 

We now ask the question of how the electron spin sys- 
tem responds to a time- varying inhomogeneous magnetic 
field. Suppose B{z,t) = {Bq + Biz)exp{—iujt). We 
show the results for diffusive dynamics, and solve Eq. 
^ with —dSo/dt as the source term. Linear- response the- 
ory for magnetic susceptibility for diffusive transport is 
well knowr£3; here we illustrate it for the specific bound- 
ary conditions of our model of spin separation, where 
various time scales discussed above will be manifest on 
the frequency domain. In response to the oscillating field, 
spin density changes as S' = [Ri{uj)+R2{z, uj)] exp(— iwi), 
where 



Ri - 

i?2 = 



Bo 
BiL 



1 



-iujTi + 1 ' 



ri>0 



cos(g„z) DqlTi + 1 



(24) 
(25) 



Spin relaxation is primarily taken over by Ri , which mea- 
sures the response of the uniform (q = 0) components of 
the spin density. On the other hand, R2 collects the terms 
responsible for diffusion, as diffusion modes {q > 0) are 
the first ones to achieve equilibrium (^t <C Ti). The total 
response, calculated for our model system, is displayed in 
Fig. ||. At small frequencies spins can adiabatically (in 
equilibrium) follow the local and instantaneous B{z,t): 
at lowtj, S « {Bq/ BiL + z/ L) exp{—iujt). At greater fre- 
quencies, first the spin-relaxation peak (in the 3i?i) and 
shoulder (in 5Ri?i) appear, while at w ~ l/ir, a second 
peak and shoulder (now due to R2) appear, as the time 
scale of the diffusive regime is reached. The second peak 
signals the SG effects, where dissipation is due to drift 



spin currents; S reaches negative values of —1/2, showing 
spin separation (compare with Fig. ^). Spectral response 
of the spin current is Js{t) = {l—dR2/dz) exp(— iwt), and 
shows a structure only around uj « I/^t, as spin current 
relaxes during the transit time (and not on Ti scale) . 

Finally, we discuss some issues related to a possible 
experimental observation of our findings. What we call 
a SG-like effect is an effective spin (not particle) separa- 
tion of electrons in metals and semiconductors. Let us 
summarize the time scales involved. Ballistic transport 
lasts for femtoseconds up to a picosecond, diffusive tran- 
sit across a micron sample can take from a picosecond to 
a nanosecond, and spin relaxation times caoJpe between 
a fraction of a nanosecond to a microsecondli3. Ordinary 
SG fails to work with electrons because transverse mag- 
netic fields (say, By = —Biy) give rise to the Lorentz 
force which makes, say, moving to the left spin up elec- 
trons turn around and move to the_ught, smearing out 
spin separationO (see, however, Ref£3). In our case the 
time scale of the Larmor precession, ti, = m/eB, is large 
enough to be neglected. Indeed, for such a large Bi as 10 
T/cm, the transverse field would be of order 10 Gauss for 
a micron sample, turning an electron around in ~ 5 
ns, long after spin separation sets in. In addition, orbital 
effects are inhibited due to momentum scattering. One 
can still go to a one- or two-dimensional saHwle to study 
SG with ballistically propagating electronai3, if Larmor 
precession is faster than momentum relaxation. 

A SG-like spin separation should be observable in both 
metals and semiconductors. (We are not aware of any 
experimental method of measuring directly the spin cur- 
rent, in our case the diffusive ptegime, although theo- 
retical proposals exist-see Ref.Ea). One way of mea- 
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suring a nonequilibrium spin in metalSpis the Silsbee- 
Johnson method of spin-charge couphngij. One can ei- 
ther switch an external inhomogeneous magnetic field, or 
inject nonequilibrium spin into a metal in a static field, 
to measure the time evolution of the spin. This can be 
accomplished, for example, by placing a ferromagnetic 
electrode on the top of a sample's edge, and measure 
the voltage across the iiiterface, which is proportional to 
the nonequilibrium spintj (the spin in the ferromagnetic 
electrode can be considered to be in equilibrium, since 
it relaxes much faster than in the nonmagnetic sample). 
The voltage would be present even when the spin cur- 
rent vanishes (that is, in the quasiequilibrium regime), 
and thus can be monitored with the sub-Ti (not tr) res- 
olution. Gradient Bi can create spin at the sample edges 
of about fisBiL/ Ep spins per electron [Ep is the Fermi 
energy), which, for typical values of, say, _Bi w 1 T/cm, 
L — 10 /im, and Ep w 10 meV gives about 1 spin per 
10^ electrons (note that L must be smaller than the spin 
diffusion length to observe the separation). For compar- 
ison, in the Johnson-Silsbee spin injection experiment 1 
spin per 10^^ electrons was detectedliil. 

In semiconductors like GaAs, the traditional tool to 
observe spin polarization of the carrier has been pho- 
toluminescence polarization detectiorO. A degenerate 
semiconductor with Ep w 1 meV (and otherwise the 
same conditions as above) would be polarized to about 
0.01% (for more sensitivity a greater Bi or L, or a ma- 
terial with a larger factor could be used), emitting 
light with circular polarization of the same order. If 
the sample is n-doped, for example, and the edges form 
the interfaces with a p-doped material, the spin polar- 
ization of the light emitted at the edges would be op- 
posite in the quasiequilibrium regime, demonstrating the 
SG-separation. Pico-to-micro-second resolved pump and 
probe photoluminescence measurements in an inhomoge- 
neous but static magnetic field could follow the evolution 
from ballistic regime to full equilibrium of a semicon- 
ductor spin system, yielding information not only about 
spin, but also about charge transport, as seen from our 
calculation. In addition to the optical technique and the 
Johnson-Silsbee method, one could also in principle ob- 
serve our predicted effect using the magnetic resonance 
force microscopy. 

In summary we have studied transient spin dynam- 
ics of itinerant electrons in metals and degenerate semi- 
conductors placed in an inhomogeneous magnetic field. 
In particular, we have solved exactly a spin dynamics 
model based on the Boltzmann equation, and demon- 
strated that the spin evolution proceeds through four 
distinct modes: ballistic, diffusive, quasiequilibrium, and 
equilibrium. An effective spin separation is possible in 
the quasiequilibrium regime, where the spin current van- 
ishes and the spin is in equilibrium only with the inho- 
mogeneous component of the magnetic field. 

We thank Igor Zutic and Xuedong Hu for useful dis- 
cussions. This work was supported by the U.S. ONR and 
DARPA. 



APPENDIX A: 

To demonstrate the effect of spin-flip scattering on 
charge and spin transport, consider electrons in a sim- 
ple metal, scattering elastically off impurities at the rate 
WkA,k'A'- Spin-flip events are characterized by Wkt,k'i 
and result mainly from the spin-orbit interaction. The 
collision integral is 



k'A' 



[Wk'A',kA/k'A'(l - hx) ~ W^kA.k'A'/kA(l - /k'A')JAl) 



The first term describes scattering from k'A' to kA, which 
increases /kA, while the second term represents reversed 
processes, which decrease /kA- Factors like /kA(l — /k'A') 
ensure that the initial state is occupied and final state 
empty, in accord with the Pauli principle. 

Let the system is driven off equilibrium by an applied 
electric field E and an inhomogeneous chemical potential 
A(5/i. The latter does not disturb the equilibrium elec- 
tronic density, but maintains an inhomogeneous nonequi- 
Ibrium spin polarization (XSfi is essentially the driving 
term for spin diffusion caused, for exarnple, by spin in- 
jection). We seek the solution to Eq. |l| with the RHS 
Eq. M in the form /kA = /kA - {dfo/deu.)(t>k\, where 
now /^ = /o(ek — XSfi). After linearization the Boltz- 
mann equation becomes 

(eE + WSn) • Vk = ^ W^kA,k'A' (<^kA - 0k' A') , (A2) 

k'A' 

where we left out the Lorentz forc e as unimportant for the 
present discussion. Equation A2 has the formal solution 



fkA 



— — eE ■ ak — AV5^ • bk. 



(A3) 



with ak and bk satisfying the integral equations: 
Vk = X! [W^kT,k'T(ak' - ak) + M/kT,k'i(ak' - ak)] ,(A4) 



and 



Vk = 5] [WkT,k'T(bk' - bk) - M^kT,k'i(bk' + bk)] .(A5) 



Vectors ak and bk have magnitudes of order £, and their 
knowledge allows to calculate the tensors of charge and 
spin conductivities: a = 2e^ ^j.(— 9/^/(9ek)vkak and 
as = 2e^ ^^(— 9/^/9ek)vkbk. The spin conductivity 
as is related to the spin diffusivity Ds as as = e^gpDs, 
since the nonequilibrium spin is dS = gpSfi. Here we use 
spin conductivity instead of the more usual spin diffusiv- 
ity only to stress the contrast with charge conductivity. 

If there is no spin-flip scattering, ak = bk, the effective 
mean free paths are the same for both currents, and the 
conductivities are equal: as — a. Spin-flip scattering, 
however, implies ak ^ bk, and so it plays different roles 
in charge and spin transport. Assume, for a moment, 
that scattering is isotropic, and energy surfaces spherical. 
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Then Eqs. A4 and A5 can be solved exactlyB by intro- 
ducing transport relaxation times r and r^: at = — '''Vk 
and bk = — TgVk, and, with 9 being the angle between 
Vk and w'^, 

- = V [M^kT,k'T(l - cosd) + M^kT,k'i(l - cos^)] ,(A6) 

T ^ — ' 



and 



r = H [^kT.k'T(l - + WkT,k'i(l + COS0)] .(A7) 



In charge transport, spin-conserving and spin-flip pro- 
cesses contribute in the same way: they are weighted by 
the well known 1 — cos 0, which suppresses contributions 
from small-angle scattering as ineffective in degrading 
charge current. Spin transport is a different story. Here 
spin-flip processes come with 1 + cos 9, and backscatter- 
ing {9 K, tt) is the least effective in degrading spin cur- 
rent, while small- angle events contribute most. There 
is illuminating physics behind this: Spin up and down 
electrons move antiparallel to each other, so if any spin 
flip is accompanied by the velocity reversal, the current 
does not change. But if the velocity stays the same, 
the effect is maximal, as if the electron sp in d oes not 
flip, but the velocity reverses. Equations A6 and A7 
are not valid in more general cases of anisotropic bands 
and inelastic anisotropic scattering, but, with some cau- 
tion, they are still pa. useful approximation (justified by 
variational analysisEj) in the interesting case of electrons 
scattering by thermal fluctuations (lattice or spin) at low 
T. Such fluctuations will allow only small-angle scatter- 
ing, suppressing contributions from spin-conserving, but 
not from spin-flip processes. Spin-flip processes are thus 
much more important for spin transport than for charge 
transport. For example, the contribution of the phonon- 
induced spin-flip scattering to l/r falls as T^, while to 
1 /ts only as (this follows from Yaleikr. theoryEJ as 
also conflrmed by a numerical calculationlla) . The con- 
tribution from the spin- conserving electron-phonon in- 
teraction falls as (the Bloch-Griineisen law), for both 
charge and spin currents. 



APPENDIX B: 




FIG. 4: Integration contour for K{q, t) of Eq. |B1| with the 
horizontal and vertical axes representing the real and imagi- 
nary part of p, respectively. The integral in Eq. Bl, which 



runs along C\, is the same as the integral along the path C2 
cutting out the branch line from p- to p+, plus the residue 
at po. 



The integral can be evaluated by a suitable contour de- 
formation in the complex pla ne, as indicated in Fig. |^. 



The original integral in Eq. Bl, which goes along Ci 



is the same as the integral over C2 plus the residue at 
Pq. The path C2 cuts away the branch line extending 
from p_ = — I/tto — iqvp to p+ = —1/Tm + iqvp, from 
the complex plane. The residue is evaluated for the pole 
at Po = + qvp cot{qvFT) present for qvpT < 7r/2 

(deflning the Riemann sheet for arctan to go from —tt/2 
to 7r/2). 

The result of the contour integration can be formally 
written as 



K{q,t) = (qvpr) e 



exp [qvpt cot (qvpT)] 



sin (qvpr) 



qvpT 



where the vertical bars denote the singular (principal) 
part of the Laurent series in terms of qvpT, of the ex- 
pression inside. An alternative formulation is 



The solution of Eq. || (or, equivalently, Eq. |T^) can be 
written with the help of Laplace transform as the integral 
in the complex plane: 



K{q,t) 



dp eP\{p + tqv+l/rmr') ^g^^ 

'00-\-icr 



where u > 0. For a degenerate system considered in the 
text, 



{{p + iqv + 1/t„i) ^) = 



1 



2iqvf 



In 



p + iqvp + 1/t„ 



p - iqvF + 1/t„ 



(B2) 



Kiq,t) = -iT/t)e~'/^"^ J2 nF^{qvFt)iqvFTr, iB4) 



n— — 00 



where functions Fn (x) are the coefficients of the Laurent 
series: 



exp[xcot(y)] = ^ i^n(a;)?;". 



n— — 00 



satisfying the recursion relation 

77 + 1 

F„(x)+<'(x) = —F^+,{x), 



(B5) 



(B6) 



7 



with the boundary conditions -^^(O) = and -F/j(0) = 
(5„._i for n < —1. In principle, all the functions Fn{x), 
n — —2, .., — oo can thus be generated from 



F-i(x) — sin(a;), 



(B7) 



which is readily obtained from Eq. p36| . 

The limiting case for ballistic transport can be ob- 
tained from Eq. B4 by letting qvpT (and thus also qvpTi) 
going to infinity. The result is 



K{q,t) ^ Ki,,u{q,t) = 



qvpt 



(B8) 



On the other hand, the diffusive limit can be obtained 
by letting qvpT to zero, in which case the vertical bar 
in Eq. B3 can be removed (the singular part equals the 
whole part since the regular part vanishes). By expand- 
ing cot{qvpT) w l/{qvpT) — {qvpT)/3 and denoting as 
D = VpT/3, one obtains 



K{q,t) « KMQ,t) = cM-t/Ti - q^Dt). 



(B9) 
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